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Abstract We propose a compact atomic clock based 
on ultracold Rb atoms that are magnetically trapped 
near the surface of an atom microchip. An interroga- 
tion scheme that combines electromagnetically-induced 
transparency (EIT) with Ramsey's method of separated 
oscillatory fields can achieve atomic shot-noise level per- 
formance of 10~^^ / y/r for 10^ atoms. The EIT signal can 
be detected with a heterodyne technique that provides 
noiseless gain; with this technique the optical phase shift 
of a 100 pW probe beam can be detected at the photon 
shot-noise level. Numerical calculations of the density 
matrix equations are used to identify realistic operat- 
ing parameters at which AC Stark shifts are eliminated. 
By considering fluctuations in these parameters, we es- 
timate that AC Stark shifts can be canceled to a level 
better than 2 x 10~^^. An overview of the apparatus is 
presented with estimates of duty cycle and power con- 
sumption. 



1 Introduction 

Recent progress in miniaturizing cold- and ultracold- 
atom systems has opened up the possibility of using ul- 
tracold matter and Bose-Einstein condensates (BECs) 
for portable applications. Inertial sensing, gravimetry, 
and timekeeping all stand to benefit from recent ad- 
vances in the field HHISHISlEliaElCQ . Some key devel- 
opments include compact ultrahigh vacuum (UHV) sys- 
tems, distributed- feedback diode lasers, and silicon fabri- 
cation techniques for both trapping with atom microchips 
and MEMS-based fiber-optic components. Several of these 
technologies were recently used to demonstrate BEC pro- 
duction in a robust, portable system having a volume of 
only 0.4 m^ [10 . 

We propose a transportable, miniature atomic clock 
based on laser-cooled ^^Rb atoms that are magnetically 
trapped near the surface of an atom microchip ^. 



Cooled, trapped atoms have significant potential in ad- 
dressing needs for high-accuracy clocks having a small 
footprint. For example, Doppler shifts are significantly 
reduced by the cold temperatures, while the spatial con- 
finement provided by the trap prevents decoherence from 
wall collisions. In addition, atom trapping can in princi- 
ple make the clock insensitive to orientation in a gravi- 
tational field and can accommodate a dynamic environ- 
ment subject to modest g-forces. Moreover, implemen- 
tation with an atom microchip enables one to consider 
its integration with other chip-based cold atom sensors. 

Currently, chip-scale atomic clocks (CSACs) based 
on hot alkali vapors represent the state-of-the-art in 
miniature atomic clocks pT|[T2l [T3l. Although instabil- 
ities down to 6xl0~^^ have been achieved, there are 
obstacles to achieving better performance from these de- 
vices. First, the coherence time of the atoms is shortened 
by collisions between the atoms and the cell walls. Al- 
though the use of buffer gases reduces dephasing, coher- 
ence 

times are still limited to only a few ms. Second, this 
use of a buffer gas introduces pressure-dependent fre- 
quency shifts that are affected by temperature changes. 
Third, continuous excitation of the atoms by the laser 
field generates AC Stark shifts and linewidth broadening 
that contribute to long-term drift [14[ ll5j. 

With the limitations of clocks based on vapor cells on 
the one hand, and the unprecedented accuracies better 
than 10~^^ in trapped, cold- atom macroscale clocks on 
the other [16, IT, cooling and trapping are advantageous 
towards enabling further progress in the development of 
miniature atomic clocks. Atom trapping is crucial for 
lengthening coherence times and hence improving clock 
stability. For example, consider a physics package mea- 
suring no more than 10 cm in any dimension. Without 
trapping, the coherence time of an atom will be lim- 
ited by the travel time of ~ 50 ms before gravity forces 
it to hit a wall, where the resulting collision will de- 
stroy the coherence. The resulting linewidth is ~10Hz, 
regardless of how cold the atoms are. Permitting the 
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clock's physical dimension to be even longer than 10 cm 
in the direction of gravity, the atomic transverse temper- 
ature still limits travel times to less than 50 ms, assuming 
typical atom temperatures from a magneto-optical trap 
(- 150 /iK for s^Rb). 

For magnetically trapped, ground state Rb atoms, 
coherence times greater than 2 s have been observed ex- 
perimentally [6l[T8]. Our target for a background pres- 
sure of 10~^^ Torr and atom temperature of less than 

1 /iK should result in coherence times of at least 1 s, 
providing an interrogation time equivalent to that of 
a 1.25 m fountain clock. The resulting Ramsey fringe 
linewidth of 0.5 Hz would be an order of magnitude nar- 
rower than what can be obtained with untrapped cold 
atoms and two orders of magnitude narrower than what 
can be obtained in a vapor cell. Such improvements more 
than compensate for the loss of atoms that accompany 
cooling and trapping. 

Due to their insensitivity to magnetic fields, the 
m = magnetic sublevels normally used as the clock 
states in alkali-atom clocks cannot be used. Instead, the 
low-field seeking, magnetically trappable \F = l^mp = 
— 1) and \F = 2,mF = +2) ground state sublevels can 
be used. In a magnetic field near 3.23 G, the first-order 
Zeeman shift of this transition cancels, helping to reduce 
the sensitivity of the clock frequency to fiuctuations and 
inhomogeneities in the trapping fields [6,18 . 

For interrogating the atoms, we consider Ramsey's 
method of separated oscillatory fields with pulsed elec- 
tromagnetically induced transparency (EIT) [19 . Imple- 
mentations of Ramsey spectroscopy using pulsed EIT 
have been studied extensively in atomic beams and va- 
por cells [20,21,22 . However, to our knowledge, no at- 
tempt has yet been made to experimentally observe this 
phenomenon in cold, trapped atoms, where it has been 
proposed as a way to eliminate AC Stark shifts in atomic 
clocks [li[22l[23]. 

This paper is divided into four sections. The first 
part is a basic overview of the proposed experimental 
approach. The second part presents theoretical calcula- 
tions of Ramsey spectroscopy with pulsed EIT. From nu- 
merical solutions to the density matrix equations, clock 
stability is estimated for realistic operating parameters. 
The third part focuses on systematic effects. We show 
that for a proper choice of laser pulse durations and 
detunings, AC Stark shifts can be eliminated. Dephas- 
ing due to magnetic field inhomogeneity and mean-field 
shifts are also reviewed. In the last part, an overview of a 
proposed experimental implementation is presented with 
focus on a compact UHV vacuum system. Estimates of 
the duty cycle and power consumption for a first gener- 
ation implementation are also presented. 

2 Experimental Overview 

A simplified energy diagram of the D2 transition in ^^Rb 
is shown in Fig. IT] The clock states |1) = |5^5'i/2,i^ = 




6.8 GHz 



Fig. 1 Simplified energy diagram of ^ Rb. The clock states 
|1) = |5'5i/2,F = l,mF = -1) and |2) = |5'5i/2,F = 
2,mF = +1) can be magnetically trapped and have an 
energy spacing equal to the ^ 6.8 GHz ground-state hy- 
perfine splitting. The excited state |3) can be either the 
|5'P3/2,i^' = 2,m,./ = 0) or the \5^P3/2,F' = l,mj.. = 0> 
sublevel. 

l,mF -1) and |2) \b'^Si/2,F = 2,mF = +1) have 
an energy splitting UJ21 equal to the ~6.8 GHz ground- 
state hyperfine splitting. These sublevels were chosen as 
clock states since they are both low-field seekers that 
can be trapped in a magnetic field minimum. A pair 
of laser beams with wavelengths near 780 nm and fre- 
quencies separated by ~6.8 GHz are used to drive two 
electric dipole transitions: a cr+-polarized probe beam 
drives the transition |1) |3) and a cr~-polarized cou- 
pling beam drives the transition |2) |3). State |3) 
is one of the magnetic sublevels in the excited 5^P3/2 
manifold. Due to selection rules, this state must have 
magnetic component tuf' = and total angular mo- 
mentum F^ = F±l = lor2. Therefore the states 
\^^P3/2.F' = l.TUF' = 0) or |52P3/2,F' = 2,mF' = 0) 
can serve as |3). Unprimed labels F and m refer to 
ground states in the 5^51/2 manifold, and primed labels 
are used for excited states in the 5^P3/2 manifold. 

All atoms must be in states |1) or |2) at the beginning 
of the operational sequence (see Fig. |2|. This can be 
achieved either by optically pumping atoms into one of 
these states or simply waiting for atoms in other states, 
which are magnetically untrappable, to be ejected from 
the trapping region. Probe and coupling laser fields are 
then pulsed for duration ti, using EIT to put the atoms 
into a steady-state, coherent superposition of states |1) 
and 1 2): 

1^) = ^i^m, (1) 
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Fig. 2 The interrogation scheme uses pulsed EIT to imple- 
ment Ramsey's method of separated oscillatory fields. The 
first pulse puts the atoms into a superposition of clock states. 
Atoms then precess in the dark for an interrogation time T, 
after which the laser difference frequency is compared to the 
clock transition frequency by measuring an optical phase shift 
inscribed on the probe beam during the second EIT pulse. 



where f2p and f2c are the Rabi frequencies of the probe 
and coupling transitions, respectively. We assume f2p = 
Qci resulting in an equal superposition (|1) — |2))/v/2. In 
this case, the initial EIT pulse is functionally equivalent 
to the initial 7r/2 pulse in Ramsey spectroscopy. 

After the laser fields are extinguished, the atoms free- 
ly precess in the dark for an interrogation time T. This 
step is equivalent to the second stage in Ramsey spec- 
troscopy, yielding fringes of width 1/(2T) Hz. The linear 
decrease in linewidth with interrogation time is one of 
the great benefits of Ramsey's method. By comparison, 
spectroscopy methods that continuously probe atoms, 
such as the coherent population trapping (OPT) inter- 
rogation done in CSACs, yield clock performances that 
improve only as the square-root of time. 

The coherent superposition state created with the 
first EIT pulse is dark with respect to a pair of laser 
beams whose difference frequency is equal to {jJ2\' This 
dark state is detected by applying a second EIT pulse of 
duration ^2 and monitoring the optical phase shift of the 
probe beam. If the laser difference frequency is equal to 
a;2i, then the atomic sample remains transparent, and 
no phase shift is inscribed by the atoms onto the probe 
beam. If the laser difference frequency does not equal 
then the atoms will not be dark and the probe beam 
will pick up a phase shift. We assume that the Rabi 
frequencies of the second pulse are the same as those 
for the first pulse, otherwise the second pulse will pump 
atoms into a different superposition state, via Eqn. [l] 
leading to a non-zero phase shift of the probe beam even 
if the laser difference frequency equals uj^i- 

For magnetically trapped Rb atoms, recent theoret- 
ical work has investigated the impact of field inhomo- 
geneity and cold collisional shifts on coherence time, fre- 
quency drift, and accuracy |24j. We work within the pa- 
rameters of references [18] and [24 and consider a sam- 
ple of Rb atoms that has been evaporatively cooled to a 
temperature of = 500 nK. Atoms are trapped near the 
surface of an atom microchip in a harmonic trap formed 
from electrical currents running through conductors on 
the chip surface. The atomic density will have a Gaussian 
spatial profile with a FWHM of v'S ln2 A:BTa/(mcj2)^ 
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Fig. 3 Schematic for optical heterodyne detection of a 
100 pW beam at the shot-noise level. The probe is split into 
strong (Pg ~ 1 mW) and weak (P^ ~ 100 pW) beams. The 
weak beam is frequency shifted with an acousto-optic mod- 
ulator (AOM) before passing through the atomic cloud. The 
weak and strong beam are recombined on the face of a pho- 
todetector. The electronic beat at the AOM frequency is de- 
modulated using a phase-shifted copy of the AOM drive. The 
proper choice of phase selects between in-phase and quadra- 
ture signals, one of which is directly proportional to the phase 
shift inscribed by the atoms onto the weak beam. 



where /c^ is Boltzmann's constant and m is the atomic 
mass. For a transverse trap frequency of ijJi = 100 Hz 
and an axial trap frequency of cOa = 10 Hz, the size of 
the atomic cloud is 25 jam in the transverse directions 
and 250 /am in the axial direction. For simplicity, we as- 
sume the atomic sample to have a uniform cylindrical 
density profile with these dimensions. 

Mode-matching the waist of the probe beam to the 
very tight transverse dimension of the atomic cloud pre- 
sents a challenge for achieving shot-noise limited detec- 
tion of the probe beam. To avoid saturation of the op- 
tical transitions (i.e. Oc^Op ^ (27r)6MHz), the probe 
power must be less than 1 nW. Such low powers cannot 
be detected at the shot-noise level using conventional 
PIN or avalanche photodiodes. To help circumvent this 
limitation, we choose the \F^ = 2, mi?/ = 0) sublevel for 
the excited state |3). Compared to the \F' = I^ttif' = 0) 
sublevel, the probe transition to this state is weaker, and 
therefore larger probe powers can be tolerated. 

We present two optical-RF double heterodyne tech- 
niques for detecting laser powers as low as 100 pW at the 
shot-noise level. In the first technique (see Fig. [3|, the 
probe laser is split into weak and strong beams of powers 
Pw and Pg, respectively. The weak beam (P^ ~ 100 pW) 
is shifted by ~80MHz in an acousto-optic modulator 
(AOM) before passing through the atomic cloud. The 
strong beam (P^ ~ 1 mW) propagates parallel to the 
weak beam, but is laterally translated so that it does 
not interact with the atoms. The two beams are recom- 
bined with a beamsplitter and the heterodyne beat at 
80 MHz is detected with a silicon photodetector. The 
beat signal has an amplitude proportional to yP^^Pg, 
with the strong beam acting as a local oscillator that pro- 
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vides noiseless gain for detection of the weak beam. At 
80 MHz, technical noise on the lasers has rolled 
off and the beat can be detected at the shot-noise limit. 
The shot noise is dominated by that of the strong beam 
and is proportional to The resulting signal-to-noise 
ratio of ^/P^Ps/ = is determined by the shot- 
noise of the weak beam. The electronic signal at 80 MHz 
is demodulated with an adjustable phase (j) that selects 
between the in-phase and quadrature components. 

A limitation of this technique is that mechanical drift 
of the optics will introduce systematic errors. The suc- 
cess of this scheme, therefore, depends crucially on the 
mechanical stability of the setup. One method for avoid- 
ing this limitation is to repeat the probe measurement 
a second time without atoms. This second measurement 
calibrates the systematic error and can be used to cor- 
rect the first measurement. The use of two measure- 
ments assumes that mechanical drift is slow compared to 
the time required to perform two probe measurements. 
For mechanical vibrations on faster timescales, a dual- 
wavelength interferometer can be used. Here, systematic 
phase shifts are measured at the same time as the atom- 
induced phase shift using a laser beam at a wavelength 
much different than that of the probe. 

The second heterodyne detection technique is based 
on phase modulation of the probe beam. Here, strong 
modulation generates a weak carrier with strong side- 
bands that provide noiseless gain for carrier detection. 
In the absence of an atom-induced optical phase shift, a 
phase modulated beam will produce no electronic beat 
signal upon photodetection. However, once the carrier 
is phase shifted, the phase modulation is partially con- 
verted into amplitude modulation that can be electron- 
ically detected. The benefit of this approach is that op- 
tical phase shifts due to mechanical instabilities affect 
all frequency components identically; this common-mode 
systematic noise will therefore be rejected upon detec- 
tion. However, the drawback to this approach is the rise 
of additional AC Stark shifts due to the multiple fre- 
quency components of the probe beam. 



3 Theoretical Calculations 



In order to estimate clock stability, we numerically solved 
the density matrix equations for the three-level system. 
The solutions are used to obtain the susceptibility, from 
which the optical phase shift inscribed by the atoms onto 
the probe beam can be calculated. The phase shift is nu- 
merically integrated as a function of time to obtain the- 
oretical Ramsey fringes as well as numerical estimates 
for signal-to-noise ratio and clock stability. 



3. 1 The Dynamical Equations of Motion 

Consider the three-level system shown in Fig. [l] in the 
presence of an electromagnetic field given by 

^ = ^ [A31 exp(-zcj3it) + A32 exp(-zcj32t) + c.c] . (2) 

The field consists of two components: a probe beam of 
amplitude A31 and frequency uo^i that is nearly resonant 
with the transition |1) |3), and a coupling beam of 
amplitude A32 and frequency a;32 that is nearly resonant 
with the transition |2) |3). Light changes the optical 
properties of the medium because it affects the evolution 
of the atomic density operator p. The evolution equation 
for p is 

= - ^[^a + H,nU P] + Rp. (3) 

where Ha is the atomic Hamiltonian in the absence of 
light and Hint = —D-S{t) is the dipole interaction 
Hamiltonian between the light field and the atom. D 
is the dipole moment operator and Rp is the relaxation 
operator. 

An electric dipole transition between clock states |1) 
and 1 2) is not allowed by selection rules. The resulting 
atomic dipole operator has the form 

|3)(1|ax3i + |3)(2|ax32+c.c., (4) 

where firnn = {m\D\n) is the dipole transition matrix el- 
ement. The Hamiltonian H = Ha ^ Hint in the rotating- 
frame approximation can be written 

H = hu;i\l){l\ + hu;2\2){2\ + hu;s\S){3\ 
- ^ (r23i|3)(l|e— ^ + r232|3)(2|e— t + ^.c.) (5) 

where f2si = AX31 • A^i/h and i732 = /l^32 • ^32/^ are the 
Rabi frequencies of the probe and coupling transitions, 
respectively. 

Evolution of the electromagnetic field is governed by 
Maxwell's equations, which reduce to 

where V is the polarization of the medium. The atomic 
dipole operator is given by Eqn. |4j so the single atom 
polarization is 

P = + M23P32 + C.C. (7) 

The polarization V of the medium is given by the rela- 
tion 

V= ^(P31+P32+C.C.) (8) 
= ^a(Ml3P31 + M23P32 + C.C.) (9) 

where is the density of atoms. The steady-state re- 
sponse of the term P31 oscillates at the frequency ujsi 
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(P31 (X exp(— cjait)) and the term P32 oscillates at the 
frequency 002,2- 

Eqns.[3]and[5]show that the coherences, or off-diagon- 
al elements of the density matrix, can be written as 
P31 = cTsi exp(-zcj3it) and P32 = (J22 exp(-icj32t), where 
the functions dij are slow functions of time as com- 
pared with the exponentials. It is convenient to factor 
out the fast time dependence from the problem and work 
with slow- varying elements of the density matrix. This is 
achieved by going to the interaction representation, i.e. 
factorizing the Hamiltonian in Eqn.|5]as H = Ho-\-V and 
introducing the interaction-representation density oper- 
ator (J by the relation a = ex.-p{iHot / h) p ex.-p{—iHot / h) . 
The density operator a obeys the equation 



dt 



a = 



h 



[Vi,a]^Ra 



where 



Vi = exp{iHot/n)V exp{-iHot/n) 



(10) 



(11) 



is the interaction Hamiltonian. Similarly, R^r is the re- 
laxation operator in the interaction representation. 

The Hamiltonian Hq should be chosen so that el- 
ements of the density matrix a are slow functions of 
time, i.e. do not contain fast time scales on the order of 
inverses of laser frequencies. A convenient choice is 



Ho = h{u;s-u;si)\l){l\ 

+ fi(c^3-c^32)|2)(2| + fic^3|3)(3| 



(12) 



and 



y = fiZi3i|l>(l| + ?i^32|2){2| 



J(r23i|3){l|e-^--* 



r232|3)(2|e- 



-iuJ32t 



+ C.C.), (13) 



where A^i = uj3i — {ujs—uji) and Z\32 = C(;32 — (CJ3— C(;2) are 
the detunings of the two lasers from the frequencies of 
the corresponding transitions with which they are nearly 
resonant. 

The slowly- varying density matrix a is related to the 
full density matrix p by 



a = 



P2ie 



Pii 



Psie 



iujsit 



P32e 



P22 

iU32t 



P23e 

P33 



■ (14) 



The interaction Hamiltonian Vj in Eqn. 11 is given by 
the expression 



v, = h 



^31 -J7i3/2' 

Z\32 -/?23/2 

-f?3l/2 -i732/2 



(15) 



The relaxation matrix is given by 

RCF = 



AlCTSS -7cC^12 -7clC^13 
-7cC^21 A2Cr33 -lc2(^23 
— 7clC^31 -7c2Cr32 -^Cr33 



(16) 



The relaxation matrix includes the total spontaneous de- 
cay rate F of the excited state; partial decay rates /31 
and /32 of the excited state to the clock states |1) and 
1 2), respectively; optical decoherences 7ci and 7c2; and 
the Raman decoherence 7c. The system is assumed to be 
open, implying that atoms are lost upon spontaneously 
decaying to states other than |1) or |2). As a result, the 
total decay rate of the excited state F = (27r)6.0MHz is 
larger than the sum of the partial decay rates (see Table 
[1]). We ignore decoherence between the clock states (i.e. 
7c = 0) and set the optical decoherences equal to half of 
their corresponding partial decay rates (i.e. 7ci = /31/2 
and 7c2 = ^32/2). 

The set of equations for the elements of the density 
matrix takes the form 



d_ 
dt 
d 
dt 



dt 



(721 = {iS - 7c)cr21 + 2 [^23^^31 - ^31^^23], 



0-31 



(^23 



^(^0 + ^) -7ci 



0-31 



-[i73i(crii - ass) + ^32^^21], 



-i(Z\o--)-7 



c2 



(^23 



(17) 



"2 [^23(^^22 - Cr33) + ^13C^2l], 

rsicrs3 + - [f2is(Tsi - c.c]. 



— (722 = ^32^^33 " ^ [^32^^23 " C.C.], 



dt 
d 

dt""'' 



-FCFSS + - [-QisCTsi + ^32^^23 - Cc], 



where the Raman detuning 5 = Asi — As2 = (^31 — 
^32) — (^2 —^1) is the difference between laser difference 
frequency and the atomic clock frequency CJ21 = 0J2—^\' 
The common- mode detuning Z\o = ^As\ +Z\32)/2 is the 
average of the two laser detunings. 

The polarization of the atomic sample is given by 



V= -(P31+P32+C.C.) 



^a[Ml3Cr3ie 



-iusit 



-^23^326"'"^^'+ C.C.]. (18) 



The relative dielectric constant of the medium is e = 
I+X7 where x is the medium susceptibility. Using Eqns.[6] 
and 18, X becomes 



2naf2isha: 



31 



2na/?23ftcr32 



(19) 



where uOp and cOc are the probe and coupling optical fre- 
quencies, respectively. 
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Probe 


Coupling 




Transition 


|1)-|3) 


|2) ^ |3) 




Partial Decay Rate (MHz) 


0.51 


1.52 




Dipole Moment (Cxm) 


7.32 X 10"^° 


1.27 X 10"^^ 




Rabi Frequency (kHz) 
Intensity (W/m^) 


150 


150 




0.245 


0.082 




Power (pW) 


128 


43 


Oh 

O 

Oh 



Table 1 Transition properties used for numerical simula- 
tions p5^. The excited magnetic sub level \F' — 2^mpi — 0) 
is used for state |3). The intensities of the two lasers are cho- 
sen so that the Rabi frequencies for the two transitions are 
the same. For obtaining laser power from intensity, a beam 
waist of 25 /xm is assumed. 



S.2 Numerical Analysis 



Eqns. 17 were solved numerically using properties for the 
relevant transitions in ^^Rb (see Table [l]). Both probe 
and coupling lasers are on resonance with their respec- 
tive transitions (i.e. S = and Aq = 0). Several factors 
must be considered when choosing appropriate values for 
the Rabi frequencies. A longer response time, as needed 
for a smaller Rabi frequencies, permits the electronic sig- 
nal processing to operate at lower bandwidths, which 
simplifies the detection circuitry. On the other hand, 
smaller Rabi frequencies result in lower optical powers 
that in turn leads to lower detection signal-to-noise ra- 
tios. Taking these considerations into account, we chose 
f2c = = (27r)150kHz for numerical estimates. As- 
suming a beam waist of 25 /im, the corresponding probe 
power is 130 pW. 

Fig. [4] illustrates how the probe and coupling beams 
modify the density matrix populations. We assume the 
atomic sample starts entirely in state |1) (i.e. pn = 
1, P22 = P33 = 0), although any choice of initial condi- 
tions will lead to the same steady-state response of equal 
clock state populations {pn = P22)' The steady-state 
population P33 = in the presence of resonant laser fields 
is the clearest indication that the sample becomes opti- 
cally transparent. Note that the excited state population 
P33 briefly becomes non-zero. Although a large fraction 
of this excited state population decays into the dark co- 
herent state, the remainder decays to ground state sub- 
levels other than the clock states. Atoms in these other 
ground states will no longer be trapped and can therefore 
escape from the atomic cloud. The result is steady-state 
populations that sum to less than 1. Assuming these un- 
trapped atoms are initially at rest, either gravity or a 
magnetic force will force them to fall out of the trap in 
less than 5 ms (ignoring collisions with other atoms). 
Although these untrapped atoms remain in the atomic 
cloud for the entire duration of the first pulse, they are 
not resonant with probe and coupling beams due to the 
magnetic field that lifts the degeneracy of magnetic sub- 
levels. 
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Fig. 4 The change in atomic populations as the initial FIT 
pulse is applied. The atoms start in the atomic eigenstate 
|1) (pii = 1)- For Qp = i7c, the steady-state populations of 
the ground states are equal (pn = ^22)- The excited state be- 
comes briefly populated (^33 > 0) before attaining its steady- 
state value of 0. 
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Fig. 5 The change in ground state coherence between the 
two clock states during the initial FIT pulse. Pumping of 
the atoms into a coherent superposition state gives rise to a 
non-zero, steady-state value of p^2i • 



Fig. [5] shows how coherence between the clock states 
is created during the initial pulse. For an equal super- 
position state, the coherence between clock states, given 
by the off-diagonal element p2i of the density matrix, 
attains a final value P2i^ equal in magnitude to the clock 
state populations pn and p22- It is this non-zero coher- 
ence that differentiates atoms in a superposition state 
from atoms in equal populations of two different atomic 
eigenstates. 

After the initial EIT pulse, atoms precess in the dark 
for time T, during which the coherence p2i picks up a 
phase shift factor exp(— ^T). After precession, the sec- 
ond EIT pulse is applied. We assume that values of the 
Rabi frequencies Qc and i?^, Raman detuning 5, and 
common-mode detuning Aq are equal to those used for 
the first pulse. Eqns. [17] are again solved numerically, 
this time using P2i^exp(— (5T) for the initial value of the 
coherence. The clock state populations obtained at the 
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Fig. 6 Optical phase shift as a function of time during the 
second EIT pulse. When the Raman detuning is on resonance 
= 0), the atomic sample is transparent and the probe beam 
picks up no phase shift from the atoms. The largest phase 
shifts occur for 8 — (n/2-\- 1/4) /T, where n is an integer. 
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Raman detuning 8 (Hz) 

Fig. 7 Integrated optical phase shift as a function of Raman 
detuning 5 shows Ramsey fringes with a FWHM of Fr — 
1/(2T), where T is the interrogation time. 



end of the first pulse are used as the initial values for the 
second pulse. 

The optical phase shift is given by A(j) — 27r/(n' — 
1)/A, where / is the length of the atomic sample in the 
propagation direction of the probe laser and the real 
component of the refractive index is given by 
n' = Re[Y^l + x(^p)]- Fig. [g] illustrates the calculated 
optical phase shift as a function of time. When the laser 
difference frequency is resonant with the clock frequency 
UO21 (i.e. = 0), the sample is transparent and there is 
no phase shift (red). The maximum phase shift (black) 
occurs for 5 = ±2ii x 1/(4T) Hz. 

The detected photocurrent is given by 
/ = 21Z\/PwPs(l)^ where IZ is the responsivity of the pho- 
todetector. To maximize signal-to- noise, the photocur- 
rent is numerically integrated up to time ^2- Ramsey 
fringes of FWHM Fr = 1/(2T) are visible in Fig. [t) 
which shows the integrated photocurrent as a function of 
Raman detuning S. Unlike fringes that are obtained from 
absorption measurements, fringes obtained from optical 



Atom Number 
Atom Temperature Ta 
Atom Cloud Length / 
Atom Cloud Waist 
Peak Atom Density Ua 
First Pulse Duration ti 
Interrogation Time T 
Fringe FWHM Fr 
Second Pulse Duration t2 
Photodiode Responsivity 71 
Weak Probe Power Pw 
Strong Probe Power Ps 



10 



10^ - 10^ 
500 nK 
250 /im 
25 /xm 
- 10^^ cm" 



13 



-10 /is 

1 s 
0.5 Hz 
-10 /is 
0.5 A/W at 780 nm 
130 pW 
1 mW 



Table 2 Assumed values for operational parameters. 



phase shifts are dispersive and therefore have a zero- 
crossing at ^ = 0. As a result, for small ^, the integrated 
photocurrent is directly proportional to the correction 
that must be applied to the laser difference frequency. 



3.3 Signal-to- Noise Ratio (SNR) and Frequency 
Instability 

The expected noise is obtained by integrating the noise 
current in over the duration ^2 of the second EIT pulse. 
This yields a noise charge Assuming the photodetec- 
tion noise is determined by the photon shot noise of the 
strong probe laser, the rms noise current is given by 



(20) 



where e is the electron charge and the circuit bandwidth 
27rZ\/ = 1/ri^ is the inverse of the characteristic response 
time tr of the circuit. The noise charge qn exhibits a 
random walk of duration yielding a mean-square noise 
charge (g^) of 



{ql) = 2Dt2^ 



(21) 



where D = aiTR is the diffusion coefficient. Eqn. [21] 
assumes that the circuit response time tr is much faster 
than t2 (i.e. tr <C ^2). The rms noise charge aq can then 
be expressed 



(22) 



Note that in the limit r^^ <C ^2, o-q depends only on ^2 
and not on the circuit bandwidth. 

The signal is obtained from the slope S of the Ramsey 
fringe near ^ = 0. The signal-to- noise ratio SNR=6'/<jg 
is numerically calculated to be ~7.5xl0^ for the exper- 
imental parameters listed in Table [2] The fractional fre- 
quency instability, as limited by noise of the optical de- 
tection, is 



^(det)/ ) ^1F^1_ 

= 1.0 X 10"^^ 



(23) 



8 



Daniel M. Farkas et al. 



where /o = cc;2i/(27r) = 6.8 GHz is the clock frequency 
and r is the integration time. The SNR of the optical 
detection is almost two orders of magnitude larger than 
the SNR due to atom shot noise alone (SNR~ = 
10^ for Na = 10^ atoms): 



(atoms) 

y 



(r) 



= 7.4 X 10"^"* 



T. 



(24) 



Aside from systematic effects and technical noise sources, 
clock performance will therefore be limited by the atom 
shot noise. 



4 Systematic Effects 

In this section, two systematic effects that impact the 
accuracy and long-term stability of the clock are con- 
sidered. The first, AC Stark shifts, can be eliminated 
with a proper choice of laser parameters (e.g. intensities, 
EIT pulse duration, and common- mode laser detuning). 
The second set of shifts arise from the combined effects 
of magnetic field inhomogeneity and mean field shifts. 
The sensitivity of the clock frequency on relevant ex- 
perimental parameters, including laser intensity, optical 
frequency, pulse timing, trap magnetic fields, and atom 
number is presented. The overall clock performance is 
calculated based on conservative estimates of drift in 
these parameters. 



4.1 AC Stark Shifts 

AC Stark shifts, or light shifts, are one of the most prob- 
lematic sources of systematic error in both miniature 
and full-size atomic clocks. The term "light shifts" refer 
to any effect that modifies the clock frequency through 
the atomic interaction with the laser fields (and other 
oscillating electromagnetic fields). The result is a clock 
whose accuracy ultimately depends on external experi- 
mental parameters instead of the inherent properties of 
the atomic species. In addition to accuracy, drift in these 
parameters can affect clock stability. 

The advantage of pulsed EIT as a clock interroga- 
tion scheme is that a proper choice of Rabi frequencies 
(as determined by the laser intensities), EIT pulse dura- 
tion, and common-mode laser detuning leads to a clock 
that is free of all light shifts. Our goal in this section is 
not to provide an accurate calculation that includes all 
sources of light shifts, which would be a computation- 
ally intensive task. Rather, we show that the shift-free 
configuration can be obtained for realistic operating pa- 
rameters by just accounting for the the largest source of 
non-resonant light shifts. However, since the identifica- 
tion of shift-free parameters is best done experimentally, 
our numerically studies strongly suggest that an experi- 
mental approach would succeed in achieving a shift-free 
clock. 



Light shifts can be divided into two categories. The 
first are shifts that arise from the resonant interaction 
of the two laser fields with the three-level atomic system 
(i.e. the idealized system represented in Fig. [T]). These 
shifts, which we label "resonant", have been studied ex- 
perimentally in atomic beams and have been found to 
arise when the initial EIT laser pulse is extinguished be- 
fore the steady-state response has been reached. For the 
case of a closed atomic system with excited state de- 
cay rate the resonant light shift AS can be expressed 
mathematically as 



-tan 



Dn 



1 



1 



sm{f3ti) 



(25) 



where Do is the difference in the initial clock state pop- 
ulations PS]. The Raman damping rate a and Raman 
dispersion f3 are given by 



1 n'^r 



2r2 



4Z\2 



r2 



(26) 
(27) 



where f2 is the Rabi frequency of both probe and cou- 
pling transitions (i.e. f2p = f2c) and the interrogation 
time T is assumed to be much longer than the pulse du- 
ration ti. The assumption of small common-mode detun- 
ings (Z\o < r) yields a ^ Q^jilF) and f5 ^ Q^A^jF^. 
The additional assumption that the argument of the arc- 
tangent is much less than 1 allows us to simplify Eqn. 25 
to 



T r2 



(28) 



e 2r 



Eqn. 28 shows how the resonant shift depends on var- 
ious experimental parameters. The shifts are inversely 
proportional to interrogation time T, reflecting the fact 
that the longer the atoms precess in the dark, the less 
the clock transition is perturbed by the lasers. The ex- 
ponential term arises from the extinction of the first EIT 
pulse before the steady-state response is reached. This 
term is reduced by increasing t\ or i?, which allows the 
system to reach equilibrium more quickly. In addition, 
shifts are reduced if the initial clock state populations 
are closer to their steady-state values. As a result, the 
shift is proportional to I^o- Last, the resonant light shift 
can be made positive or negative via the linear depen- 
dence on Z\o. 

Additional AC Stark shifts arise from the presence of 
other atomic states. We label these shifts "non-resonant" 
since the lasers that couple these transitions are detuned 
by many linewidths. In general, these shifts scale in- 
versely with laser detuning, so the dominant contribu- 
tion will come from other nearby magnetic sublevels in 
the 5^6'i/2 ground state and 5^P3/2 excited state hy- 
perfine manifolds (see Fig. [s]). Although there are 24 
sublevels in total, only three of the eight sublevels in the 
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Fig. 8 Energy diagram of the 5^ 5*1/2 ground state and 
5^i^3/2 excited state hyperfine manifolds in ^^Rb. Resonant 
(solid line) and non-resonant (dashed line) transitions al- 
lowed by the probe (coupling) laser are shown in red (blue). 
Only transitions that couple to the clock states in the ground 
state manifold are considered. 



ground state can be magnetically trapped; contributions 
from the five other untrappable ground-state sublevels 
are ignored. Furthermore, the \F = 2,mF = 2) ground 
state sublevel, which can be magnetically trapped, can- 
not be populated via spontaneous decay of an = 
excited level, leaving the two clock states as the only 
ground state sublevels that need consideration. Of the 
sixteen excited state sublevels, only six can be non-reson- 
antly coupled to the clock states with circularly polarized 
light. 

The excited non-resonant sublevel that produces the 
largest AC Stark shift is |4) = = 1, m^' = 0). There 
are two reasons for this. First, the state |4) is the only 
non-resonant excited sublevel that couples to both laser 
fields. Second, non-resonant transitions to excited sub- 
levels with \m'\ = 2 are coupled with a laser detun- 
ing of ^ 6.8 GHz, while the detunings for transitions to 

= sublevels are given by the smaller 100 to 300 MHz 
splittings of the excited 5'^Ps/2 manifold. In fact, for our 
choice of state |3), the detuning of ~ 157 MHz for state 
1 4) is smaller than that of all other = non-resonant 
sublevels. 

To cancel light shifts in the clock, a resonant light 
shift that is equal in magnitude, but opposite in sign, 
to the sum of all non-resonant light shifts is intention- 
ally introduced; the resulting resonant and non-resonant 
light shifts sum to zero. To demonstrate this cancelation 
numerically, the density matrix in Eqn.[T4]was enlarged 
to a 4x4 matrix for the inclusion of state |4) in the sys- 
tem dynamics. The expanded density matrix equations 
were solved and the integrated optical phase shift was 
calculated as a function of 5 over the center part of the 
central Ramsey fringe. The clock frequency is given by 
the value of S at which the fringe signal equals 0, which 
was found numerically for different parameters. 
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Fig. 9 The total light shift as a function of common-mode 
detuning Aq for different values of the pulse duration ti. 



Several assumptions were made in solving the 4-level 
system. First, the slaved-state approximation was made 
by setting the time-derivative of the coherence between 
the two excited states equal to 0. This is equivalent to 
assuming that the population of state |4) is negligible. 
Second, the resonant light shifts in Eqn. [28] depend only 
on the duration ti of the first FIT pulse. Non-resonant 
light shifts, on the other hand, arise for any time the laser 
fields are present, and are therefore affected by the dura- 
tion of both pulses. To simplify the calculations, the first 
and second pulse durations are assumed to be equal, or 
ti = ^2- Third, the Rabi frequencies are assumed to both 
equal a fixed value f2c = (27r)150kHz, and the ini- 

tial atomic population is assumed to be entirely in state 
|1) so that 1^0 = 1. As a result, the resonant light shift 



in Eqn. [28] is modified by changing the pulse durations 
and common-mode detuning Aq. Fourth, both resonant 
and non-resonant light shifts scale inversely with inter- 
rogation time T in the limit T ^ ti] the value T = 1 s 
was used for all calculations. 

The dependence of the clock frequency on pulse du- 
ration is presented in Fig. [9] which shows the total light 
shift scaling linearly with Aq. For values of ti = t2 be- 
tween 4 and 12 /is, the total light shift can be eliminated 
for detunings Aq between -4 and -14 kHz. The steepest 
slope of -0.34 mHz/ (kHz T) was found for ti ~ 6 /is. 
This worst-case dependency of the clock performance on 
fluctuations in Aq corresponds to a fractional frequency 
shift of -5.0 X 10-^^. 

The common- mode detuning Af at which light shifts 
cancel is plotted as a function of pulse duration in Fig.[TO 
The maximum value of Af occurs near 6 /is. At this 
point, the clock will be least sensitive to fluctuations 
in pulse duration. However, this point also corresponds 
to the greatest sensitivity to fluctuations in Aq. There- 
fore, pulse durations different from 6 /is should also be 
considered, even though this will increase sensitivity to 
timing jitter. 

Clock shifts arising from fluctuations in the laser 
parameters were estimated by numerically calculating 
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Fig. 10 Shift-free value of the common- mode detuning 
as a function of pulse duration ti. 



the total light shift for small changes of these parame- 
ters. Table [3] summarizes the results for ti = 6 /is and 
ti = 20 /is. The units are normalized to typical uncer- 
tainties of the corresponding parameters (i.e. a 1 kHz 
uncertainty in Z\o, a 0.1% shift in ti, etc.). Even though 
we assume the first and second pulses are equal in du- 
ration, timing jitter affects light shifts differently for the 
two pulses. Therefore, we consider the effects of tim- 
ing jitter for each pulse separately. Changes in the Rabi 
frequencies arise from laser intensity variations and are 
expected to be the largest source of laser-induced in- 
accuracy for ti = 20 /is. For ti = 6 /is, the dominant 
sensitivity is due to Z\o, whose fluctuations depend on 
how stable the optical frequencies are. For both values 
of ti, the clock is least sensitive to timing jitter in the 
laser pulses. The overall performance is slightly better 
for ti = 20 /is due to the reduced sensitivity to Aq. In 
this case, all inaccuracies due to laser-induced effects are 
less than 2 x 10~^^. Active stabilization of the laser in- 
tensity could reduce this even further. 

The dependency of the clock frequency on the Rabi 
frequencies, even under light-shift-free conditions, im- 
plies that the Af cannot be experimentally identified 
by varying probe and coupling laser intensities. Instead, 
Af can be found by measuring the clock frequency as a 
function of Aq for different values of T. Since both res- 
onant and non-resonant light shifts scale inversely with 
T, the clock frequency varies linearly with Aq with a 
slope proportional to 1/T. When plotting the shifts for 
all values of T on the same graph, the lines correspond- 
ing to different values of T will intersect at the shift-free 
detuning Af (see Fig. 3c in Reference [19]). 



4.2 Dephasing and Coherence Times 

After the atoms have been pumped into a coherent su- 
perposition state, two effects lead to dephasing and sys- 
tematic frequency shifts: differential Zeeman shifts and 
mean- field shifts. These two effects generate a spread in 
clock frequencies Au that cause the atoms to precess at 



Parameter 


ti = 6 /is 


ti 20 /is 


units 


Ao 


-5.0 X 10-^^ 


-7.0 X 10-^^ 


(kHz T)-i 


ti 


+5.9 X 10"^^ 


-1.7 X 10"^^ 


(0.1% T)-i 


t2 


-6.0 X 10"^'^ 


-4.2 X 10"^^ 


(0.1% T)-i 




+2.6 X 10-^^ 


-8.2 X 10-^^ 


(0.5% T)-i 




+5.4 X 10"^^ 


+1.3 X 10"^^ 


(0.5% T)-i 



Table 3 Light shifts arising from uncertainties in the 
common- mode detuning Z\o, pulse durations ti and t2, and 
Rabi frequencies i7c and i7p. Values are obtained for ti = 6 /iS 
and ti = 20 /is assuming i7c = = (27r)150kHz, ti = t2, 
and Aq = Af. Shifts are inversely proportional to T. Units 
are expressed in terms of typical uncertainties of the corre- 
sponding parameters. 



different rates, washing out the Ramsey fringes. These 
two effects can be made to largely cancel each other for 
a proper choice of trap magnetic fields. 

Differential Zeeman shifts arise from the inhomoge- 
neous magnetic field of the trap. In a harmonic trap, hot- 
ter atoms with greater kinetic energy are able to travel 
farther from the trap center, where they are exposed 
to larger magnetic fields and therefore greater Zeeman 
shifts. The magnetic field of the trap can be written 



B = Brr 



at {x^ 



y 



(29) 



where Bm is the minimum value of the field at the center 
of the trap; x, and z are the trap's spatial coordinates; 
and at and ai are the field curvatures in the transverse 
and longitudinal directions, respectively. In the presence 
of a weak magnetic field, the potential energy U (x, z) 
of an atom scales linearly with the field: 



U (x, y, z) = qf rnp Mb [Bn 



+ at [2 



y 



- aiz 



(30) 



where and tuf are the ^-factor and magnetic quan- 
tum number of the atomic state, /i^ is the Bohr magne- 
ton, and m is the mass of the atom. The product gr^F 
is equal to 1/2 for both clock states. 

The Breit-Rabi formula predicts that the linear Zee- 
man shifts of the two clock states are equal at Bq ^ 
3.23 G. Near the second-order Zeeman shift of the 
clock transition can be expressed 



Au, = b{B - B^f 
where h = 431Hz/G^ and Vrn 



(31) 

4497 Hz. Plugging 
Eqn. [30] into Eqn. [31] gives the position-dependent Zee- 
man shift Avz{x^ y^ z) of the clock transition 



Avz{x,y,z 
b[B, 



Bo + at{x^ + y^) + aiz^] ^ - (32) 



The Zeeman shift varies as the fourth power of position. 
More importantly, the sign of Bm — Bq can be changed 
to make the curvature at the trap center positive {Bm > 
Bo), negative {Bm < Bq), or flat {Bm = Bq). 
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Parameter 


Sensitivity 


Parameter Uncertainty 


Frequency Uncertainty (mHz) 


Fractional Frequency Uncertainty 




-35 Hz/G 


10-^ = 30 /xG 


-1.1 


1.5x10-^^ 


Ta 


1 Hz//iK 


10-^ = 5pK 


0.005 


7.4x10"^^ 


N 


-5 /xHz/atom 


Vn = 10^ 


-5 


7.4x10"^^ 


f 


0.3 mHz/% 


0.5% 


0.15 


2.2x10"^^ 



Table 4 Differential Zeeman shifts and mean- field shifts arising from uncertainties in the trap bias field Bm, atom temperature 
Ta, atom number and population partition /. Typical uncertainties in the parameters, and the corresponding shifts in the 
clock frequency are listed. Sensitivities in the second column are taken from reference [24j . 



For cold, non-condensed atoms in the s-wave regime, 
atom- atom collisions generate a frequency shift that de- 
pends on scattering length and atomic density. For the 
clock transition, this mean-field shift can be expressed 

^ ^ ^ 

= Q^12<^12(^1 - ^2) + a22<^22^2 " <^11<^11^1 

m 



= — n[a22<^22 - <^ii<^ii 



m 



/(2ai2ai2 - aiiaii - a22a22)], 



(33) 



where m is the atomic mass, ni and 712 are the atomic 
densities for states |1) and |2), n = ni + 712 is the to- 
tal atomic density, and the partition / = (ni — n2)/n is 
the fractional change in population densities due to un- 
equal superpositions. We assume the atoms are pumped 
into equal superposition states (i.e. / = 0). The s-wave 
scattering lengths are represented by a^j and the zero- 
distance two-atom correlation functions are given by a^j. 
For non-condensed atoms, an = 0^22 = <^i2 = 2. Us- 
ing the measured scattering lengths of a22 = 95.47ao, 
ail = 100.44ao, and ai2 = 98.09ao pF, where ao is the 



Bohr radius, Eqn. 33 can be simplified to 



Avc{x,y,z) = -4.97^aon(x,?/,2;), 
m 



(34) 



which shows that the shift is both negative and propor- 
tional to atom density n{x^y^z). Due to the Gaussian 
density profile of the ensemble, atoms at the center of 
the trap experience a more negative shift than atoms at 
the edge of the trap. 

The mean-field shift can be mostly canceled by the 
differential Zeeman shift when is chosen to be slightly 
less than Bq. In this case, the Zeeman shift Au^ has pos- 
itive curvature at the trap center, in contrast to the neg- 
ative curvature of the collisional shift. Taken together, 
the combined shift is nearly flat across the spatial extent 
of the trap, thus reducing the spread of clock frequencies 
and increasing coherence times. 

A theoretical framework has been developed in ref- 
erence [24 for better estimating dephasing times of Rb 
atoms in a magnetic trap. In this work, the total shift 
Avz^ was either spatially averaged over a Maxwell- 
Boltzmann distribution of atoms of finite temperature or 
averaged over the motion of the atoms in the trap. The 
resulting average shift was used to compute the variance 
cr^, from which the dephasing time is given by (V^ttct)"^. 



Calculations showed good agreement with previous ex- 
perimental data plITS]. 

This framework was used to both minimize the line 
spread of the clock transition and to increase the de- 
phasing time by varying atom temperature, bias mag- 
netic field, peak atom density, atom number, population 
partition, and trap frequencies [24 . Dephasing times as 
long as 20 s, corresponding to a line spread of 11 mHz, 
were obtained for experimental parameters similar to 
those used in references [6 and [18 : a temperature Ta = 
500 nK, peak densities of ~ 2.0 x 10^^/cm^, a total atom 
number N = 10^, a population partition / = 0, axial 
trap frequencies between 1 Hz and 10 Hz, and transverse 
trap frequencies between 300 Hz and 1 kHz. 

The effects of fluctuations in magnetic field, tem- 
perature, atom number, and partition on clock perfor- 
mance are summarized in Table |4j with values for the 
sensitivities coming from reference [24 . The largest frac- 
tional frequency sensitivity is due to fluctuations in atom 
number N . Unlike drift in the other parameters, fluc- 
tuations in atom number have a white noise spectrum, 
and therefore this source of clock instability will improve 
as 1/Vt for repeated clock cycles (assuming each clock 
cycle starts with the same nominal number of atoms). 
Although this systematic effect will slightly degrade sta- 
bility, it will not contribute to the noise floor that will 
ultimately limit the long-term stability of the clock. 

The stability of the trap bias field Bm is limited by 
the stability of the electrical current flowing through 
the atom chip; a 10~^ relative stability can be readily 
achieved. Since the final atom temperature is determined 
by the RF evaporation frequency and bias field Bm^ the 
temperature should have the same relative stability as 
Bm- One source of drift in the population partition / is 
fluctuations in the Rabi frequencies of the two optical 
transitions. For probe and coupling laser intensities sta- 
ble to 1%, the Rabi frequencies, and therefore the parti- 
tion /, will be stable to 0.5%. The effects of fluctuations 
in the trap frequencies are negligible [24J. 

Given that the magnetic trap is expected to be the 
largest source of systematic error, we briefly consider a 
miniature clock that confines atoms optically (e.g. an 
optical dipole trap) . This approach could require several 
Watts of laser power, however this is not necessarily a 
larger power drain than what is already needed for mag- 
netically trapping on a microchip. Similar to a fountain 
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clock, this configuration could use the m = ground 
state sublevels for the clock transition, helping to reduce 
sensitivity to stray magnetic fields. Due to the lack of a 
magic wavelength for alkali atomic clocks based on op- 
tical dipole traps, an additional light shift that depends 
on the intensity and wavelength of the trapping laser 
will arise [27 t [28 | l29l . Although further calculations for 
the case of Rb are needed, this light shift could in prin- 
ciple be eliminated either using the pulsed EIT scheme 
presented here [19] or with the addition of an external 
magnetic field [ 291130( 131] . The result would be an atomic 
clock free of both light shifts and Zeeman shifts. 



5 Description of Proposed Apparatus 

In this section, we present a overview of how recent tech- 
nological developments in ultracold atom technology can 
be adopted for the creation of a miniature, transportable 
atomic clock. There are four critical elements to the pro- 
posed system: the miniature ultrahigh-vacuum (UHV) 
cell with integrated atom microchip; magnets; the laser 
system, optics, and detection system; and control elec- 
tronics with a low phase noise microwave oscillator. Our 
focus here is primarily on the UHV cell and laser system. 
The optical detection system has already been discussed 
above (see Fig.[3|. Miniaturization of control electronics 
and low phase noise microwave oscillators has already 
been developed for CSACs, and our proposal relies on 
the same technology. 

The clock will be based on a double MOT vacuum 



cell system (see Fig. 11) that has been used for several 
cold-atom experiments and applications in our lab at the 
University of Colorado at Boulder. The double MOT 
configuration offers a substantial advantage in provid- 
ing high atom fiux, fast loading times, and long MOT 
and magnetic trap lifetimes. With an internal volume 
of 30 cm^, this cell has a vacuum level better than 
5 X 10~^^ Torr and has been used to create Bose-Einstein 
condensates near the surface of an integrated atom mi- 
crochip with a duty cycle as fast as 3 s [10,. The mi- 
crochip is anodically bonded onto the end of one of the 
chambers, where it forms a wall of the cell. Current is 
passed into the vacuum through UHV-compatible elec- 
trical vias fabricated on the chip. Low vacuum levels 
are maintained by a 21/s ion pump and non-evaporable 
getters. Cell construction utilizes silicon-to-glass anodic 
bonding and standard metal, glass, and glass-to-metal 
construction methods. These construction techniques en- 
able high-temperature bake-outs up to 300° C, which is 
essential for good UHV performance. 

The cell consists of two chambers joined by a 0.75 mm 
aperture that permits differential pumping between the 
two chambers. In the lower chamber, a 2D+ MOT pro- 
duces a beam of cold atoms, which feeds the upper cell, 
where the atoms are captured in a 3D MOT. The 2D+ 
MOT produces a captured atom fiux of 0.5 to 1 billion 



atoms/s with less than 120 mW of total optical power. 
The 2D+ MOT chamber is operated near 10-'^ Torr of 
rubidium pressure for high capture efficiency, while the 
3D MOT chamber is maintained at a much lower pres- 
sure for long trap lifetimes. Typical trap lifetimes are 10 
to 20 s for the 3D MOT, and 4 to 7 s for a loffe-Pritchard 
magnetic trap. 

After capture in the 3D MOT, the magnetic fields 
are turned off and a few ms of polarization gradient 
cooling (PGC) is applied to bring the atom tempera- 
ture below 50 /iK. The atoms are then optically pumped 
into a single, trappable magnetic sublevel. External mag- 
netic field coils are used to transport atoms a distance of 
~ 1 cm in 300 to 500 ms, leaving the cold atomic cloud 
within 0.5 mm of the vacuum side of the atom microchip. 
Atoms are transferred onto the chip by simultaneously 
ramping up current fiowing on the chip and ramping 
down current fiowing through the external field coils. 
The chip conductor forms a "Z" configuration that, in 
conjunction with a pair of external bias fields parallel to 
the chip surface, creates a loffe-Pritchard trap approxi- 
mately 100 /im from the chip surface. Evaporative cool- 
ing with an RF antenna reduces the atom temperature 
below 1 /iK. 

There are several benefits to using atom microchips. 
First, the significantly higher trap frequencies that can 
be obtained (1-2 kHz in the transverse direction) com- 
pared to macroscopic magnetic traps speeds up rether- 
malization during evaporative cooling, helping to shorten 
the duration of this stage. Second, less RF power can be 
used during evaporation by placing the RF antenna on 
the back-side of the chip, where it is less than 0.5 mm 
away from the atoms. In addition, the ability to gener- 
ate very stable electrical currents may lead to trap fields 
that are more stable than what can be obtained with 
permanent magnets, which have at best a temperature 
coefiicient of 0.03%/°C. 

One drawback to using an atom microchip is power 
consumption. The chip trap must be on during evapora- 
tive cooling and clock interrogation, for a total of several 
seconds. Assuming a wire resistance of 2 1] and a steady- 
state current of 1.5 A, the resulting power dissipation is 
more than 4 W (ignoring the additional power needed to 
generate external bias fields). One way to mitigate this 
loss is to use a planar configuration of permanent mag- 
nets on the chip as the final trap [32 . Chip wires can 
then be used for the much shorter stage of transferring 
atoms to the chip. 

Distributed feedback (DFB) diode lasers can be used 
to create a compact laser system that, compared to ex- 
ternal cavity diode lasers, is significantly more robust 
towards mechanical vibrations. In addition to the one or 
two lasers needed for EIT, up to two additional lasers 
are needed for the MOT: one for the \^Si/2^F = 2) ^ 
\5Ps/2^F' = 3) cooling transition and the other for the 
|5S'i/27^ = 1) ^ 1^^3/2^ = 2) repumping transition. 
MOTs of 1 to 2 billion atoms have been achieved in the 
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Fig. 11 Schematic of the double MOT vacuum cell system. 



two-chamber cell of Fig. 11 using 120 mW of combined 
cooling and repump power. However, to obtain 1 million 
atoms in a magnetic trap after evaporative cooling, we 
estimate only 100 to 200 million atoms will be needed in 
the 3D MOT. As a result, we expect only 20 to 30 mW 
of laser light will be necessary. By tuning the power and 
frequencies of the cooling and repump lasers, these lasers 
could also serve as the FIT probe and coupling lasers. 
This option would reduce volume and power consump- 
tion by minimizing the number of lasers in the system. 



Stage 


Duration (ms) 


Power 


MOT CooUng 


1000 


0.1 


PGC 


5-10 


0.1 


Transfer to Chip 


300-500 


5 


Evaporative Cooling 


1000-4000 


6 


First FIT Pulse 


-0.01 


5 


Interrogation 


1000 


5 


Second FIT Pulse 


-0.01 


5 



Table 5 Timing and power consumption estimates for one 
operational cycle of the proposed clock. 



6 Timing and Power Consumption 

Table[5]lists timing and power consumption estimates for 
one operational cycle of the clock. The total duration 
is between 3 and 6 s, and the clock instabilities calcu- 



lated in Section 3.3 must be increased by the square-root 
of this time. The dominant power drains in the physics 
package are the RF signal used for evaporative cooling 

1 W); the power dissipated by the chip wires 5 W) 
during evaporative cooling and clock interrogation; and 
the power used to drive external magnetic field coils. 

Operation of the cold- atom clock in conjunction with 
a CSAC provides a means of duty-cycled operation to re- 
duce overall power consumption. For example, the cold- 
atom apparatus can be operated once every hundred 
seconds, providing periodic frequency corrections to a 
continuously-running CSAC. Since the power duty cycle 
is reduced to 1%, the average power consumed by the 
cold-atom physics package is reduced by two orders of 
magnitude. 



7 Conclusion and Outlook 

We have investigated the feasibility of a transportable, 
miniature atom clock based on laser cooling and mag- 
netic trapping of ^^Rb atoms near the surface of an atom 
microchip. A Ramsey sequence based on pulsed FIT, 
in conjunction with a heterodyne optical detection ap- 
proach that can detect optical powers as low as 100 pW 
at the shot noise level, can achieve detection almost two 
orders of magnitude more sensitive than performance 
limited by the shot-noise of 10^ atoms. By construct- 
ing the clock around a compact, UHV double MOT cell, 
interrogation times up to 1 s should be attainable, result- 
ing in Ramsey fringes as narrow as 0.5 Hz. Through the 
use of an integrated atom microchip, this compact ge- 
ometry both improves clock stability by increasing duty 
cycle and reduces power consumption by minimizing the 
volume of external magnetic field coils. 
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As a clock interrogation scheme, pulsed EIT permits 
implementation of a novel light-shift cancelation scheme. 
Here, resonant light shifts arising from the laser cou- 
pling of the clock states to a common excited state are 
intentionally introduced in order to cancel non-resonant 
light shifts arising from all other atomic states. The den- 
sity matrix equations were solved with the inclusion of 
a fourth atomic state that generates the largest non- 
resonant light shift; cancelation of light shifts was numer- 
ically demonstrated for realistic values of the EIT pulse 
duration and common- mode laser detuning. Estimating 
the effects of fluctuations in laser intensities, pulse tim- 
ings, and common-mode detuning on clock performance, 
residual light shifts below 2 x 10~^^ should be achievable 
assuming reasonable values for drift and uncertainty in 
these parameters. 

Clock stability will be limited by shot-noise of the 
atoms, which affect both the SNR and the combined ef- 
fects of differential Zeeman and mean field shifts. Addi- 
tionally, a slight reduction in stability arises from a duty 
cycle less than 100%. Taking all these effects into ac- 
count, a fractional frequency instability of 1 x 10~^^/v^ 
should be achievable. The accuracy of the clock should 
be limited by differential Zeeman and mean field shifts 
to better than 1 x 10~^^. 
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